Precision Measurement of the Neutron Beta-Decay Asymmetry 
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A new measurement of the neutron /3-decay asymmetry Ao has been carried out by the UCNA 
collaboration using polarized ultracold neutrons (UCN) from the solid deuterium UCN source at 
the Los Alamos Neutron Science Center (LANSCE). Improvements in the experiment have led to 
reductions in both statistical and systematic uncertainties leading to Ao — — 0.11954(55) s tat. (98) sys t., 
corresponding to the ratio of axial-vector to vector coupling A = Pa/sv = —1.2756(30). 



Precision measurements of neutron /3-decay are an es- 
sential ingredient in understanding the electro-weak in- 
teraction in the light quark sector. In particular the 
axial- vector weak coupling constant, gA, is an important 
input to understanding the spin and flavor structure of 
the nucleon [TJ [5] and is being actively studied in de- 
tailed lattice QCD calculations [3j 0]. It also plays an 
important role in a variety of astrophysical processes, in- 
cluding solar fusion cross sections important for energy 
and neutrino production in the sun [S]- 

The angular distribution of emitted electrons from de- 
cays of a polarized neutron ensemble can be expressed 
as [B] 



W(E) oc 1 + -(P)A(E) cos9, 
c 



(1) 



where A(E) specifies the decay asymmetry for electron 
energy E, v = j3c is the electron velocity, (P) is the 
mean neutron polarization, and 9 is the angle between the 
neutron spin and the electron momentum. The leading 
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order value of A(E), Aq, can be expressed as 
-2(A 2 -|A|) 



A = 



1 + 3A 2 ' 



(2) 



where A = gA/gv is the ratio of the vector to axial- vector 
weak coupling constants. Combining gA with indepen- 
dent measurements of the Fermi coupling constant Gp, 
the CKM matrix element V u d, and the neutron lifetime 
r n allows a precision test of the consistency of measured 
neutron /3-decay observables [Jj. 

The UCNA (Ultra-Cold Neutron Asymmetry) exper- 
iment is the first experiment to use ultracold neutrons 
(UCN) in a precision measurement of neutron decay cor- 
relations. Following the publication of our earlier results 
([IHS])j the UCNA collaboration implemented a num- 
ber of experimental improvements that led to reductions 
in both statistical and systematic uncertainties. These 
improvements, described below, include enhanced UCN 
storage, improved electron energy reconstruction, and 
continuous monitoring of the magnetic field in the spec- 
trometer. 

The UCNA experiment ran in 2010 using the "thin 
window geometry D" as described in [71 IS] , and collected 
a total of 20.6 x 10 6 /3-decay events after all cuts are ap- 
plied. We use the UCN source [ID] in Area B of the LAN- 
SCE. UCN are polarized by a 6 T pre-polarizer magnet 
and a 7 T primary polarizer, coupled to an adiabatic fast 
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passage (AFP) spin flipper to control the spin state [TT] . 
Upstream of the pre-polarizer magnet, a gate valve sepa- 
rates the UCN source from the experimental apparatus. 

Polarized UCN enter the superconducting spectrom- 
eter (SCS) [12], and are confined in a 3 m long, 12.4 
cm diameter diamond-like carbon (DLC) coated Cu tube 
(decay trap) with 0.7/xm thick mylar endcaps. The in- 
side surface of each endcap is coated with 200 nm of 
Be to contain the neutrons. A 0.96 T magnetic field 
is oriented parallel to the decay trap, along which de- 
cay electrons spiral toward one of two identical electron 
detector packages. Between the decay trap and the de- 
tectors, the magnetic field expands out to 0.6 T, which 
reduces the electrons' transverse momenta and pitch an- 
gles, decreasing backscattering from the detectors. Each 
detector package consists of a 16 x 16 cm 2 low-pressure 
multiwire proportional chamber (MWPC) 13J backed by 
a 15 cm diameter plastic scintillator, whose scintillation 
light is detected by four photomultiplier tubes (PMTs). 
Each MWPC has 6 ^m mylar windows at the front and 
back that separate the chamber gas (100 torr neopen- 
tane) from the spectrometer vacuum and PMT housing 
(< 100 torr N2). Cosmic-ray muon backgrounds are iden- 
tified by a combination of plastic scintillator veto paddles 
and sealed Ar/ethane drift tube assemblies [32] around 
the electron detectors. 

A typical run unit consists of a background run (gate 
valve closed), a /3-decay run (gate valve open), and a 
UCN depolarization run (see below). To partially cancel 
drifts in background and detector efficiency, we alternate 
the order of the /3-decay and background runs, and or- 
ganize the asymmetry measurements into octets with a 

spin flipper on (+), off (— ) sequence of H 1 h H — 

or hH 1 h, chosen randomly. 

Scintillator event triggers are formed by requiring at 
least 2 out of 4 PMT signals over threshold in either of 
the scintillator detectors. Due to the low mass of the 
MWPC, requiring an off-line coincidence between the 
MWPC and the scintillator rejects > 99% of external 
7-ray background. Energy deposition in the MWPC is 
calibrated against our Monte Carlo simulation to aid in 
classification of backscattering events. Cosmic-ray muon 
backgrounds are measured and vetoed off-line by requir- 
ing coincidences between any of the muon detector com- 
ponents and the electron detectors. 

Electron positions at the MWPC are determined to 
<2 mm based on the distribution of charge on two per- 
pendicular cathode grids in the MWPC [T3]- A fiducial 
cut of r < 50 mm (projected to the 0.96 T decay trap re- 
gion) is placed on the trigger side to reduce background 
and to eliminate electrons that could strike the decay 
trap walls. 

The equilibrium UCN polarization that develops dur- 
ing each /3-decay run is measured using the spin flipper to 
selectively unload polarized and depolarized UCN from 
the decay volume immediately following the /3-decay run. 
This is accomplished in two stages: first, the guide serv- 
ing as input to the 7 T primary polarizing field is switched 



TABLE I. Polarizations obtained from the two data sets: 
2010/ , which includes all 2010 depolarization and reload runs 
prior to the pump failure along with all 2009 reload runs, 
and 2010//, which includes all depolarization and reload runs 
obtained after the pump failure. 



Data Set 


(P) Polarization 


2010/ nipper off 


1.001(2) stat (5) sys 


2010/ flipper on 


0.990(l) stat (5) sys 


2010// flipper off 


0.992(5) s tat(8) sys 


2010// flippper on 


0.988(4) stat (3) sys 



to guide neutrons towards a 3 He UCN "switcher detec- 
tor" [T5] ~0.75 m below the beamline while the gate valve 
is simultaneously closed and proton pulses are discontin- 
ued. This cleaning phase, which lasts 25 s, produces a 
signal in the UCN detector proportional to the number of 
correctly polarized UCN present in the experimental ge- 
ometry at the end of the /3-decay measurement interval. 
Following the cleaning phase, the state of the spin flipper 
is changed, preventing any remaining correctly polarized 
UCN in the decay trap from exiting the geometry and al- 
lowing incorrectly polarized UCN remaining downstream 
of the spin flipper to pass through the 7 T polarizing field 
and be counted. Counting during this unloading phase 
is performed for ^200 s in order to measure incorrectly 
polarized UCN as well as background. The primary sys- 
tematic uncertainty in these measurements comes from 
any remaining correctly polarized UCN upstream of the 
spin flipper at the moment its state is changed; these 
UCN are not prevented from reaching the UCN detector 
during the unloading phase and produce a background 
whose size is of the same order as the incorrectly po- 
larized signal. Correction for this reloaded population is 
accomplished using ex situ measurements ( "reload" mea- 
surements) in which the spin flipper state is toggled for 
3 s during the middle of the cleaning phase in order to 
selectively enhance the signal from the reloaded popula- 
tion. The measured polarization in the case of a spin- 
flipper-off /3-decay run also requires correction for spin 
flipper inefficiency, which is determined using the differ- 
ence between polarizations observed for spin-flipper-off 
and spin-flipper-on along with Monte Carlo calculated 
scaling factors. Further small corrections for UCN pop- 
ulations detected in the switcher detector with low ef- 
ficiency are estimated via Monte Carlo. Based on the 
global agreement between Monte Carlo simulations and 
data, an uncertainty of 30% is attributed to all polariza- 
tion Monte Carlo calculations. An analysis of our fitting 
procedure to the switcher detector signal during depolar- 
ization runs also contributes to the systematic error. This 
includes sensitivity to the fitting intervals, along with the 
internal consistency of the extracted time constants. 

Midway through the 2010 run, a vacuum pump failure 
unexpectedly vented the spectrometer, producing pin- 
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hole leaks in the MWPC windows. For a brief period 
of operation before the windows were replaced, neopen- 
tane leaking from the wirechambers into the vacuum may 
have permanently contaminated the UCN guide surfaces, 
resulting in a change to the UCN transport characteris- 
tics of the system (e.g. a 35% reduction of UCN storage 
lifetime in the decay trap was observed after the pump 
failure). Since this incident potentially altered the equi- 
librium UCN polarization in the decay volume, separate 
polarization analyses for the periods before and after the 
pump failure were required. In order to improve the 
statistics, and because there were no observable changes 
to the experimental geometry between the 2009 run cy- 
cle and the pump failure, the set of reloaded population 
measurements obtained in 2009 was combined with the 
2010/ data acquired prior to the pump failure. The polar- 
izations determined from the "before" and "after" data 
sets are shown in Table HI 

Reconstructed event energies E Tecon are measured us- 
ing the signals from the four PMTs attached by light 
guides to the scintillator disk in each detector. The 
position dependence of light transport to each PMT is 
mapped out by filling the spectrometer volume with 
neutron-activated Xenon. Natural isotopic abundance 
Xe gas is let into the volume normally containing the solid 
deuterium UCN source, and irradiated for a few minutes 
with the source flux of spallation neutrons to produce 
a variety of radioactive Xe isotopes by neutron capture. 
After pumping the activated Xe out of the source vol- 
ume, controlled amounts are introduced into the spec- 
trometer volume. By observing the decay spectrum fea- 
tures (mainly the 915 keV /3-decay endpoint from 135 Xe 
J n = | + ) as a function of position using the MWPC, 
the position-dependent light transport of the beta scin- 
tillators is mapped out. The increased statistics avail- 
able from the Xe data compared to the previous method 
of mapping position dependence using neutron /3-decay 
data allows for increased resolution and decreased statis- 
tical noise in the position-dependent response. 

The energy response and linearity of each PMT is cal- 
ibrated with conversion electron sources ( 139 Ce, 113 Sn, 
and 207 Bi) inserted into the center of the decay trap. 
Energy losses due to the sealing foils of each source were 
determined using a collimated 241 Am alpha source and a 
silicon detector. Energy losses to 5485.6 keV alpha parti- 
cles passing through the mylar sealed source foils indicate 
an effective thickness of 9.5 /urn (compared to the nomi- 
nal 6 /im thickness specified by the manufacturer, likely 
due to the adhesives sealing the source package), uniform 
to <2% over position on the foil. 

Since the data taken for the previous publication [9] , 
which required correction for nonlinearity in some of the 
PMTs' response due to damage from sparking in PMT 
bases run at sub- atmospheric pressures, the bases as well 
as the PMTs have been replaced. The new PMTs (Hama- 
matsu R7725) show a linear response at the level of < 1%. 

The improved linearity and reduced uncertainty in po- 
sition response and source foil energy losses allow an over- 



all reduction in energy reconstruction systematic uncer- 
tainty to approximately half of the previous limit 19] . 

Variation in PMT and electronics gain is continuously 
monitored with a newly installed 207 Bi "pulser" source, 
based on the concept of [16], consisting of a scintillator 
block containing a small amount of 207 Bi mounted on the 
face of each PMT alongside the light guide from the main 
scintillator disk. A high-threshold single-PMT trigger 
distinguishes pulser events from beta scintillator events 
which typically distribute light between several PMTs. 
The ~lMeV 207 Bi conversion electron line provides a 
consistent peak for tracking gain changes, with sufficient 
statistics to measure each PMT's gain to < 0.3% over five 
minutes. Over longer time scales, the 207 Bi pulser signal 
peak was observed to drift on the order of 1% per week 
relative to periodic calibrations with the conversion elec- 
tron sources. Longer term gain stabilization for /3-decay 
data is implemented by fixing the neutron /3-decay spec- 
trum endpoint averaged over each octet (~ 8 hours) of 
runs to the expected (Monte Carlo) value, while using 
the 207 Bi pulsers to monitor and correct shorter time- 
scale drifts. 

The majority of the /3-decay events are single detector 
triggers. However, due to electron backscattering, ~ 3% 
of the events trigger both scintillators, while ~ 2.5% are 
detected by both MWPCs but trigger only one of the 
scintillators. In the first case, the initial direction of the 
electron can be determined by the relative timing of the 
triggers, while in the second case a cut based on the en- 
ergy loss in the trigger side scintillator and MWPC yields 
an identification efficiency of ^80% based on a Monte 
Carlo simulation. 

In our previous publications [7HS] the uniformity of the 
magnetic field in the decay region was checked with an 
NMR probe translated through the field with the cen- 
tral decay trap removed. Thus the field could only be 
measured at the beginning and end of a long data-taking 
period, leading to additional systematic uncertainty due 
to possible variations between the measurements. In the 
present data run the field is continuously monitored by an 
array of sixteen Hall effect sensors placed just outside of 
the decay trap, allowing the field uniformity to both be 
optimized when ramping the spectrometer magnet and 
continuously monitored. Monte Carlo simulations of the 
impact of field non-uniformity on the measured asymme- 
try provide a correction similar to analytical estimates of 
magnetic mirroring for high pitch angle events. 

In addition to the ambient backgrounds (measured 
with the UCN gate-valve closed) which are subtracted 
run-by-run, neutron captures in the vicinity of the de- 
tectors can create prompt 7's and delayed /3-decay elec- 
trons, generating an irreducible background in the exper- 
iment. We deduced a neutron generated background rate 
of ~ 0.025 Hz (under the ~ 25 Hz /3-decay rate) via excess 
events above the neutron decay endpoint. Comparison of 
observed excess event energy and position distributions 
with Monte Carlo simulations indicate the excess is con- 
sistent with a combination of UCN capture on the alu- 
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FIG. 1. Top: background subtracted electron energy spec- 
trum, combining both detector sides and spin states, over- 
laid with the Monte Carlo prediction. The measured back- 
ground spectrum is also shown. Middle: Aq vs. Brecon, shown 
with statistical error bars, and fit to a constant from 220 to 
670 keV. Bottom: corrections and their uncertainties (band) 
excluding polarization and theory contributions; positive sign 
indicating a larger |j4q|- 



minum surfaces of the detector and on the scintillator 
disk resulting from a small fraction of UCN that escape 
from a monitor hole in the decay volume as well as small 
gaps in the UCN guide system. This excess contributes 
a correction and uncertainty to the measured asymmetry 
of +0.01(2)%. 

For each run, events are binned based on reconstructed 
energy (10 keV bins) and initial direction. The rates in 
the two detectors are then computed based on the ex- 
periment live time. We applied separate spin-dependent 
blinding factors to the two detector rates, effectively 
adding an unknown scaling factor to the measured asym- 
metry that was constrained to be within 1.00(5). After 
determination of all cuts, corrections, and uncertainties, 
this factor was removed. For each /3-decay/background 
run pair, the background rate is subtracted from the j3- 
decay-run rate bin by bin. The reconstructed energy 
spectrum (background subtracted, averaged over the two 
spin states) is shown in Panel (a) of Fig. [T] overlaid with 
the measured background (SignaLBackground ~ 124 be- 
tween 220 and 670 keV). Also overlaid is the Monte Carlo 
predicted reconstructed energy spectrum, with detec- 
tor response effects (energy resolution, trigger efficiency, 
etc.) taken into account. 

In each measurement unit (octet), a ratio of count 
rates is constructed, leading to a "super-ratio" (as de- 
fined in [5]), from which the asymmetry is determined. 



The final measured asymmetry is the statistical combi- 
nation of all asymmetry sub-units therein. 

To extract Aq, we first divide the raw measured 
asymmetry by |/3 in each energy bin to remove the 
strongest energy dependence. As described in [S] , two 
scattering related effects dominate subsequent system- 
atic corrections: the residual backscattering correction 
Abackscattcring and the angle effect A ang i c . In addition to a 
small correction due to incorrect identification of the ini- 
tial electron direction for the measured electron backscat- 
ters (where both detectors observe the electron), there 
are corrections for backscattering from the decay trap 
windows and the front windows of the MWPC that can- 
not be identified experimentally. Angle effects arise from 
the fact that the energy loss of an electron in the thin 
windows is strongly angle dependent. Low-energy, large 
pitch angle electrons are more likely to fall below the 
scintillator threshold, leading to a suppression of the ac- 
ceptance at large angles. Both of these effects were evalu- 
ated with two independent Monte Carlo simulation pack- 
ages: Penelope [17] and Geant4 [18] (version 4.9.5, us- 
ing the "Livermore" low-energy EM physics model 19i ) . 
The two simulations were benchmarked against the mea- 
sured backscattering distributions for the different types 
of backscattering events using both neutron /3-decay elec- 
trons and conversion-electron sources. The resulting cor- 
rections are shown in Table [TT] For all analysis choices 
(inclusion/exclusion of backscattering event types), the 
correction calculated from the two Monte Carlos agreed 
to within 15%. Based on observed differences between 
the simulations and the detectable backscattering data 
(e.g. two scintillator triggers and two MWPC hits for 
single scintillator triggers), we assign a fractional uncer- 
tainty of 25% to the backscattering and angle effect cor- 
rections. 

Additional theoretical contributions (beyond the sim- 
ple v/c term) must be incorporated in order to convert 
the observable neutron beta decay asymmetry A(E) to 
the underlying parameter Aq. Recoil-order contributions 
to A(E) were calculated within the context of the Stan- 
dard Model according to the formalism of [2"0lI23] . and 
the radiative correction contribution was calculated ac- 
cording to [2U2S] Q 

Applying all corrections mentioned above, the ex- 
tracted A is plotted against E Iecon in Panel (b) of Fig[l] 
Energy-dependent corrections (backscattering and angle 
effects) and their uncertainty are indicated in the figure. 



The estimated radiative correction in 1241 . Eq. 15, is based on an 
energy-independent analysis that integrates total counts across 
the whole spectrum. The "Fermi function" weighting of the spec- 
trum towards lower energies (and lower asymmetry), represented 
by the Coulomb terms 2n 2 /3~ 1 in 1241 Eq. 14, dominates the 
correction. For an analysis that extracts Aq as a function of 
energy, the bin-by-bin energy-dependent correction has the op- 
posite sign. Our previous Aq measurement [9j did not account 
for this. Updating the result with the value from table|TT]modifics 



the result from [5] to A = -0.11942 ± 0.00089 



+0.00H 
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TABLE II. Summary of corrections and uncertainties as % 
of Aq. "+" corrections increase \Aq\ from the observed un- 
corrected value. 



Systematic 


corr. (%) 


unc. (%) 


Polarization 


+0.67 


±0.56 


^backscattering 


+1.36 


±0.34 


Wangle 


-1.21 


±0.30 


Energy reconstruction 




±0.31 


Gain fluctuation 




±0.18 


Field non-uniformity 


+0.06 


±0.10 


CMWPC 


+0.12 


±0.08 


Muon veto efliciency 




±0.03 


UCN-induced background 


+0.01 


±0.02 


^statistics 




±0.46 


Theory contributions 


Recoil order |20H23] 


-1.71 


±0.03 


Radiative I21I231 


-0.10 


±0.05 



The final Aq is obtained from an average over an energy 
range of 220 to 670 keV, which was chosen, before un- 
winding the asymmetries, in order to minimize the com- 
bined statistical and systematic uncertainties. In the 220 
to 670 keV range, fitting the 10 keV binned values of Aq 
to a constant value yields x 2 /ndf = 41.7/44 (based on 
statistical error bars). The energy-averaged A is also 
very stable for different energy ranges, remaining con- 
stant within ±0.15% for ranges out to 100 to 800 keV 
(where x 2 / n df = 68.2/69). 

The uncertainties and systematic corrections to Aq are 
summarized in Table [TTJ The measured result is Aq — 
— 0.11954(55) s tat. (98) sy st. where the first uncertainty is 
statistical and the second systematic. Based on Eq. ([2]), 
we can also determine A = gA/gv = —1-2756(30). The 
present result is shown in Fig. [2] compared with previous 
high precision (a a/ A < 2%) results. 

In summary we have measured the polarized neutron 
decay asymmetry with UCN resulting in a fractional pre- 
cision of < 1%. When combined with our previous pre- 
cision result [9] with the updated radiative contribution, 
we obtain a UCNA value of A = -0.11952(110) and 
A = —1.2755(30). The consistency of our results with 
the most recent measurements from the Perkeo collabo- 
ration [30] [31], which have significantly smaller correc- 
tions compared to the pre-2000 results, may suggest that 
the uncertainties were under-estimated in some of these 
earlier experiments. This consistency of the most recent 
values of A in the context of light quark decay parameters 
is shown in Fig. [3] 

With considerable efforts underway world-wide to im- 
prove the precision of angular correlations measurements 
sensitive to lambda using cold neutron beams [3"3"rJ37| . 
there remains significant motivation to continue efforts 
to further refine corresponding measurements with UCN. 
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FIG. 2. Ideogram of values for Ao from this work (filled 
square) and recent measurements (open circles) [71 l9l I27H31] . 
The shaded band indicates the current PDG 2012 average 
value of Ao = —0.1176(11) [32], which includes the results 
of [Tll§l l27H30| . but does not include [!5T] or the work reported 
here. To account for correlated systematic errors in sequential 
measurements, the ideogram (solid curve) was constructed us- 
ing the combined result from [30] and |31) of —0.11951(50) 
reported in [31] , and the combined result of [3 [9] and this 
work of —0.11956(110), as discussed in the text. 
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FIG. 3. The light quark weak coupling Kid vs. A. Kid = 
0.97425(22) from 0+ -> 0+ decays and the neutron lifetime 
r n = 880.1(1.1) s are from PDG 2012 [32]. Values of A are the 
UCNA result from this paper, and the Perkeo II combined 
result A = -1.2755(13) from [3J. 
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